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The quantum resonances occurring with 5— kicked parti- 
cles are studied with the help of a fictitious classical limit, 
establishing a direct correspondence between the nearly reso- 
nant quantum motion and the classical resonances of a related 
system. A scaling law which characterizes the structure of the 
resonant peaks is derived, and numerically demonstrated. 
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Resonances are a widespread phenomenon both in 
quantum and in classical physics. However, quantum 
resonances are often purely quantal phenomena, not di- 
rectly related to classical nonlinear resonances. Atom 
optics has made crucial aspects of the classical-quantum 
correspondence accessible to direct experimental observa- 
tion. An outstanding example is the paradigmatic Kicked 
Rotor (KR) Model, which served as a prototype of quan- 
tum and classical Hamiltonian chaos over almost three 
decades. This model has been experimentally realized by 
a technique first introduced by Raizen and co-workers [1] , 
allowing for observation of its main physical properties. 
These are connected with the long-time evolution, and 
depend on the arithmetic nature of the effective Planck's 
constant (which is the Planck constant divided by a char- 
acteristic action) [2] . If the latter is sufficiently irrational, 
the rotor's energy saturates in time. This effect is known 
as dynamical localization. If the effective Planck's con- 
stant is a rational multiple of 47r (= Anr/q, with r, q 
integers) unbounded growth of energy typically occurs. 
These are known as quantum resonances and have no 
counterpart in the corresponding classical system. Be- 
sides confirming theoretically known facts, experiments 
have also introduced variants of the basic model, which 
have enriched the theory. A striking example was the 
experimental discovery of unexpected "quantum acceler- 
ator modes" due to effects of gravity [3] . 

Even in cases when no major variants are introduced, 
experiments inevitably involve deviations from the ide- 
alized theoretical model. The most obvious of these is 
that atoms move in lines and not in circles; though to 
some extent trivial, this fact imposes certain modifica- 
tions on the standard KR theory, especially in the case of 
the quantum resonances [4]. Another experimental lim- 
itation is that experiments are confined to not too long 
times. Consequently, they can neither detect the high- 



order (high— g) quantum resonances, nor address the ba- 
sic, as yet unsolved, mathematical problem of what de- 
gree of irrationality is required for localization. Never- 
theless, for this reason they have encouraged fruitful the- 
oretical analysis of the intermediate time regime, which 
had rarely been investigated in fine detail. 

In this letter we present experimentally relevant results 
concerning the manifestations of quantum resonances at 
finite observation times. Our analysis enlightens the 
connections between quantum and classical resonances. 
When the energy of the kicked atoms is measured after 
a fixed time, and the result is plotted versus the kicking 
period, resonances manifest in the form of peaks at the 
resonant values of the period (see Fig.l). We describe 
the line shape of such resonant peaks, and its depen- 
dence on the observation time and on the other relevant 
parameters. Our main result is a single-parameter scal- 
ing law for such peaks, which in particular implies an 
inverse square dependence of the peak width on time, 
and is probably related to the "sub-Fourier resonances" 
recently detected in a different atom-optical context [5]. 
Our derivation is based on a quasi-classical analysis of 
the purely quantal resonances and their vicinity. This 
seemingly self-contradictory task is accomplished by es- 
tablishing a direct correspondence between the quantum 
resonances and the classical resonances of a related sys- 
tem, which is not obtained in the usual classical limit of 
vanishing Planck's constant, but rather in the limit when 
the effective Planck's constant approaches its resonant 
value, which is of order unity. In this limit, the detuning 
of the driving period from resonance plays the role of the 
Planck constant. The same technique was used in [4] to 
explain the quantum accelerator modes observed in the 
presence of gravity. Here it is applied in the gravity-free 
case, leading to a different physical scenario. It may also 
be adapted to the important case when decoherence is 
added [6]. 

The specific system we consider is described, in dimen- 
sionless units, by the Hamiltonian 

= ^p^ + k cos{x) ^ 6{t'~t), (1) 

t — — OG 

where x is the coordinate, p is its conjugate momentum, t' 
is the continuous time variable, while t is an integer that 
counts the number of kicks. If x is an angle variable, con- 
fined to the interval [0, 27r), then (1) is the Hamiltonian 
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of the KR and p is the angular momentum; when quan- 
tized, it takes integer values only. The two dimensionless 
parameters that control the dynamics arc the effective 
Planck's constant r and the strength of the kick k. The 
first (second) is proportional (inversely proportional) to 
Planck's constant h, therefore kr = K is independent of 
h. K is the stochasticity parameter of the Standard Map, 
which rules the classical dynamics defined by (1). The 
energy levels th? /2 [n integer) of the free rotor (fc = 0), 
and differences between them, are integer multiples of 
r/2, while the frequencies of the driving potential are in- 
teger multiples of 27r, consequently quantum resonances 
of the kicked rotor (1) are found if t/Att is rational [7]. 
These resonances correspond to the Talbot effect in optics 
and 47r corresponds to the Talbot time [8]. For most ra- 
tional t/A-k they result in ballistic growth oip (quadratic 
growth of energy) with time. 
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FIG. 1. (a) mean energy after t = 30 kicks, vs. the kicking 
period r, for an ensemble of 10'' 5— kicked atoms, witli Gaus- 
sian initial momentum distribution (a ~ 2.7) and k = O.Stt. 
(b) quantum data are compared with the mean energies of 
an ensemble of 10® particles (circles), evolving under the 
e— classical dynamics (5). The small peak on the right of the 
resonant spike is marked by an arrow for reference to Fig. 2. 
(The peaks appear lower in (a) because the used computar- 
tional grid in t had a lower resolution than in (b).) 

If X is not a coordinate on a circle but on the full line 
instead, then p is linear momentum and (1) yields the 
Kicked Particle (KP) model which is the one relevant for 
recent atom optics experiments. Since the driving po- 
tential is periodic in x only transitions between values 
of p that differ by an integer are allowed. Consequently, 
the quasimomentum (QM) (3 (the fractional part of p) is 
conserved [4,6]. Standard Bloch theory allows for reduc- 
tion of the KP dynamics to the dynamics of a bundle of 
rotors. At any given (and constant in time) value of /3, 
the dynamics induced by (1) is formally that of a rotor 
and the one-period Floquet operator for this rotor is [4]: 
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where 9 = .Tmod(27r), and Af = —\d/d6 is the angular 
momentum operator. If t/At: is rational, the KP dy- 
namics may exhibit asymptotic quadratic growth of the 
energy only for finitely many special values of the QM 
(3. The particle is then in a spatially extended state [6]. 
Wave-packet dynamics, or incoherent ensemble averag- 
ing over momentum as is the case in experimental re- 
alizations [9] enforces linear growth of the energy as a 
function of time for rational t /An [4.6]. Away from reso- 
nances, the growth of energy in time is bounded owing to 
dynamical localization [2], so a scan of the particle energy 
vs. r after a fixed evolution time exhibits peaks at the 
resonant values of r, as shown in Fig. 1. Such peaks rise 
and shrink as the observation time is increased. Their 
height is proportional to time, and their width to inverse 
square time, as our present analysis will show. 
We restrict to the main resonances and to their vicinity, 
that is to r = 27r€ -|- e where ^ is a positive integer and e is 
small compared to 27r. We denote k = k\e\, and / = |e|7V. 



Then, using the identity e"^^" 
rewritten as [4,6] 



eq. (2) may be 



^ p-ifecos(e")/|e| ^-mp/\e\ 



where 



Ul3{t) = e 



Hpil, t) = -sign(e)/2 + /(7r£ + r/3) . 



(3) 



(4) 



If ]e| is regarded as Planck's constant [10], then (3) is the 
formal quantization of either of the following classical 
maps: 



It+i = h + ksm{et+i) , 

et+i = et±It+Tre + T0 mod(27r) 



(5) 



where ± is the sign of e. The small— |e| asymptotics 
of the quantum /?— rotor is thus equivalent to a quasi- 
classical approximation based on the e— classical dynam- 
ics (5). Here the term "classical" is not related to the 
h ^ limit but to the limit e ^ instead, so the term 
"e— classical" will be used in the following. Note that 
the effective Planck's constant r is not sm,all, and impor- 
tant quantum symmetries are taken into account in the 
e-classical description. The efficiency of this approxima- 
tion is demonstrated in Fig. 1(b), where results of exact 
quantum simulations are compared to results obtained by 
(5), for an initial incoherent ensemble of particles, each 
in a momentum eigenstate (a plane wave) and with a 
Gaussian distribution of the momentum p = nQ+ (3. For 
each particle in the ensemble, the map (5) with [3 equal 
to the QM of the particle was used to compute a set of 
trajectories started at 7 = no|e| with homogeneously dis- 
tributed 9 e [0, 27r). The final energies e-2/2/2 at t = 30 
of the individual trajectories were averaged over 9, (3, no 
with the weights imposed by the initial ensemble. The 
average energy (Et) = e~'^{If) /2 is plotted vs. t = 27r + e 
in Fig. 1(b), along with results of quantal computations 
of the ensemble averaged energy. 
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We shall show that the resonant peaks are determined 
by e— classical phase space structures. It is convenient 
to change variables to J = ±/ + 7r€ + t/3, = + 
7r(l — sgn(e))/2, thus turning the maps (5) into a sin- 
gle e-classical Standard Map (eSM), independent of the 
value of p: 



Jt+i = Jt + fcsin('!?t+i) , = 'dt + Jt 



(6) 



It will turn out that: (i) The resonant value r = 
27r£ corresponds to integrablc e— classical dynamics, and 
the resonant values of the QM /3 correspond to the 
e— classical resonant values of the action J; (ii) for small 
|e|, the e— classical dynamics is quasi-intcgrable, and the 
growth of the quantum energy is dominated by the main 
e— classical resonant island around J = 27r; (iii) at any 
time t, the ratio between the energy and its value at e = 
is a scaling func tion , notably it is a function of the single 
variable x = t^/k\e\ and not of the variables t, k, e sepa- 
rately. 

We assume for simplicity an initially flat distribution of 
p G [0, 1]; then Iq = 0, and Jq = iri + tI3q with /3o uni- 
formly distributed in [0, 1). Without loss of generality 
we also consider i = 1. Hence if |e| <^ 1 then Jq is uni- 
formly distributed over one period (in action) (tt, Stt) of 
the eSM. The resonant QM value is /3 = 1/2 [6]. Since 
Jt = ±It + n + t(3, and /q = 0, the mean energy of the 
rotor at time t is: 

{Et,,) = e-^{lf)/2 = e-\{5Jtf)/2 , SJ^ = J* - Jo • 

The exact quantum resonance e = corresponds to the 
integrable limit of the eSM, where SJt = 0. However, 
(i?t,e) is scaled by e~^, so in order to compute it at e = 
one has to compute 6Jt at first order in e. This leads to: 



dJt = \e\k^sm{eo + Jos) +r{e,t) 



(7) 



s=0 



where r{e,t) = o(e) as e — > at fixed t. The particle 
energy at time t is found from (7) by taking squares, 
averaging over Oq, Jq, dividing by 2|ep, and finally letting 
e^O: 



/r2 

o) = — 



sm\Jot/2) _ k\ 
° sin2(Jo/2) ~ 4 ' • 



(8) 



The small contribution of the initial QM in the atom's 
energy was neglected. Apart from that, (8) is equal to the 
exact quantum mechanical result [6]. The integral over 
Jo in (8) collects contributions from all the invariant tori 
Jo =const. of the eSM at e = 0, but it is essentially 
determined by a small interval ~ 2-K/t of actions around 
Jq = 27r, the main e— classical resonant torus. As that 
torus is formed of (period 1) fixed points, its own con- 
tribution is quadratic in time, so the linear growth (8) 
follows. Noting that Jo = 2tt corresponds to the reso- 
nant QM: /?o = 1/2, we see that the e— quasi-classical 
approximation explains the quantum resonances in terms 



of classical resonances of the Standard Map [11]. We 
now estimate {Et,() for |e| > 0. The |e| > dynamics is 
maximally distorted with respect to the e = one for Jq 
in the vicinity of 2n7r, that is, in the very region which 
is mostly responsible for the linear growth of energy at 
e = 0. Being formed of period-1 fixed points, the Jo = 27r, 
e = invariant torus breaks at |e| > as described by the 
Poincare-Birkhoff theorem [11], giving rise to the "main 
resonance" of the eSM, which is located astride J = 2tt 
with a size SJ^es ~ 4(A;|e|)^/^ [11]. The approximation 
(7) fails quickly therein so its contribution (£'t,e)r-es in 
the mean energy has to be estimated differently. In the 
remaining part of the e— classical phase space the motion 
mostly follows KAM invariant curves, slightly deformed 
with respect to the e = ones, still with the same rota- 
tion angles. The contribution of such invariant curves in 
the mean energy is therefore roughly similar to that con- 
sidered in the integral (8). On such grounds, in order to 
estimate {Et.e) we remove from the integral (8) the con- 
tribution of the resonant action interval near Jq = 27r, 
and we replace it by: 



where 



87r 



/2 sin2(J72) 



-SJres 



(9) 



(10) 



and J' is the deviation from the resonant value 27r. 
{Et^f_)res may be estimated by means of the pendulum ap- 
proximation [11]. Near the eSM resonance, the motion is 

described (in continuous time) by the pendulum Hamilto- 
nian in the coordinates J',!}: Hres = \{J'Y + \e\kcos('&). 
A characteristic time scale for the motion in the reso- 
nant zone is tres = (fc|e|)~^^^7 which is the period of the 
small pendulum oscillations divided by 27r. One may al- 
together remove |e| from the Hamilton equations, by scal- 
ing momentum and time by factors (fc|e|)~^/^ = 4:/5Jres, 
(/c|e|)^/^ — 1/tres respectively. Therefore, 



4f) 



k\e\G{t. 



(11) 



for an ensemble of orbits started inside the resonant zone, 
where G{.) is a parameter-free function. This function re- 
sults of averaging over nonlinear pendulum motions with 
a continuum of different periods, so it saturates to a con- 
stant value when the argument is larger than w 1. The 
contribution to the total energy is then obtained on mul- 
tiplying (11) by |e|^^JJres/(47r), because only a fraction 
^ (5Jres/(27r) of the initial ensemble is trapped in the 
resonant zone. As a result 



{Et,t)r 



-^/\'"^G{t^\) 



(12) 



When 5Jres is small, sin ( J'/2) may be replaced by J /4 
in the integrand in (10), leading to 
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FIG. 2. Demonstrating the scaling (13) at t > 2it. 
Open circles correspond to different values of the parameters 
e, fc, t, randomly and independently generated in the ranges 
\ <t < 200, 0.001 < e < 0.1,0.1 < fc < 50 with the con- 
straint 0.01 < fee < 0.2. In each case an ensemble of 2 x 10® 
e— classical rotors was used to numerically compute the scaled 
energy R{t,k,e) (13), with a uniform distributions of initial 
momenta in [0,1] and of initial in [0, 27r). Full squares 
present quantum data for k = 0.8n, t = 50 and t = 200. The 
solid line is the scaling function H{x) of (13) obtained by 
computation of the functions <I>o(.'r) (dashed) and G(x). The 
dot-dashed line has slope —1. The arrow marks the value of 
the scaled detuning x corresponding to the arrow in Fig. 1(b). 



The map (5) is easily adapted to the model in the pres- 
ence of decoherence due to random momentum jumps 
induced by external noise, e.g. the Spontaneous Emis- 
sions effects used in experiments [9,12]. A scaling law is 
again valid for the resonant peaks, in the two variables 
t/tres and tc/tres where tc is the time scale associated 
with noise [6]. 

To summarize, we have exploited a correspondence be- 
tween the dynamics of a Kicked Particle near the quan- 
tum resonances r = 2Tr£ and the classical dynamics of a 
quasi-integrable system to analyze the structure of the 
experimentally observable quantum resonant peaks. For 
these we have derived a scaling law where the scaling 
variable is the ratio between the observation time and the 
characteristic time scale of elliptic motion inside the e— 
classical island. This law provides significant new infor- 
mation about the shape and the parameter dependence 
of the peaks. 
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Replacing in (9), we obtain: 
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X = tv^fcR = t/Ues ■ (13) 
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